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MEMORANDUM  FOR  Commander,  U.S.  Army  Engineer  District,  Charleston, 

ATTN:  CESAC-TS-DT  (Mr.  Wayne  Bieganousky),  69A  Hagood  Avenue,  Charleston,  SC 
29403-5107 

SUBJECT:  Seismic  Stability  of  St.  Stephen  Hydropower  Plant,  South  Carolina 


1.  Reference: 

a.  Federal  Emergency  Management  Agency  (FEMA)  Document  356,  November  2000, 
Prestandard  and  Commentary  for  the  Seismic  Rehabilitation  of  Buildings. 

b.  ERl  1 10-2-1806,  31  July  1995,  Earthquake  Design  &  Evaluation  of  Civil  Works  Projects. 

c.  Memorandum,  CECW-EG,  5  May  2000,  subject:  Seismic  Stability  of  St.  Stephen 
Hydropower  Plant,  South  Carolina. 

d.  Memorandum,  CENWD-MT-ES,  29  October  2001,  subject:  Corps  of  Engineers  Hazard 
Reduction  Program  -  Powerhouse  Evaluation  Criteria. 

2.  At  the  request  of  the  Charleston  District,  the  Engineer  Research  and  Development  Center 
(ERDC)  conducted  a  seismic  stability  evaluation  of  the  St.  Stephen  Hydropower  Plant,  SC. 
Enclosure  1  summarizes  the  details  of  this  seismic  engineering  evaluation  of  the  powerhouse 
superstructure.  A  performance-based  approach,  as  outlined  in  the  FEMA  document  356  (ref  a), 
was  used  for  the  seismic  evaluation  to  determine  the  risk  of  the  superstructure  collapse  during  a 
major  earthquake.  Collapse  prevention  performance  requires  that  collapse  of  the  structure  be 
prevented  regardless  of  the  level  of  damage  inflicted  by  the  earthquake  on  the  superstructure. 

The  superstructure  was  evaluated  for  collapse  potential  for  two  site-specific  ground  motions:  a 
2,475-year  (return  period)  probabilistic  earthquake  event  and  a  deterministic  Maximum  Credible 
Earthquake  event.  Selection  of  these  events  in  a  collapse  potential  evaluation  is  consistent  with 
FEMA  356  seismic  analysis  procedures  for  existing  buildings.  It  was  determined  that  the  St. 
Stephen  Powerhouse  superstructure  is  vulnerable  to  collapse  for  either  of  these  major  earthquake 
events.  The  superstructure  does  not  meet  the  basic  safety  objective  (of  collapse  prevention)  of 
FEMA  356.  Additionally,  it  was  also  determined  that  the  superstructure  is  potentially  vulnerable 
to  collapse  for  earthquake  hazards  with  return  periods  greater  than  500-years. 

3.  Using  the  criteria  and  evaluation  procedures  given  in  ERl  1 10-2-1806  (ref  b),  the  St.  Stephen 
Hydropower  Plant  was  deemed  by  HQ  USACE  Civil  Works  Directorate  to  not  be  a  ‘critical’ 
structure  (See  Enel  2  (ref  c).). 
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4.  The  seismic  safety  of  all  US  ACE  buildings  has  been  evaluated  under  the  National  Earthquake 
Hazard  Reduction  Program  (NEHRP),  which  is  designed  to  insure  the  safety  of  personnel  in 
USACE  and  other  Federal  buildings  in  accordance  with  PL  101-610  and  Executive  Order  12941. 
Enclosure  3  (ref  d)  is  a  draft  memorandum  from  the  NEHRP  Director,  Ms.  Helen  Peterson, 
CENWD,  to  Dr.  Ebeling,  CEERD,  stating  that  the  procedures  outlined  in  FEMA  356  are 
applicable  for  use  in  the  seismic  stability  evaluation  of  powerhouses;  the  seismic  stability 
evaluation  conducted  by  CEERD  of  the  St.  Stephen  Hydropower  Plant  is  consistent  with  the 
FEMA  356  procedures;  and  that  the  use  of  a  2,475-year  probabilistic  earthquake  event  or  a 
deterministic  Maximum  Credible  Earthquake  event  in  this  collapse  potential  evaluation  is 
consistent  with  FEMA  356  (collapse  prevention)  seismic  analysis  procedures  for  existing 
powerhouse  superstructures. 

5.  The  next  stage  in  the  evaluation  would  be  to  investigate  modifications  necessary  to  remediate 
the  St.  Stephen  powerhouse  superstructure  so  as  to  eliminate  the  potential  for  collapse  of  the 
superstructure  for  a  2,475-year  probabilistic  earthquake  event  or  a  deterministic  Maximum 
Credible  Earthquake  event. 

6.  POC  is  Dr.  Ebeling  (601-634-3458). 

3  Ends  Robert  M.  Ebeling,  PhD 

Civil  Engineer  Supervisory  Research  Civil  Engineer 

Ralph  W.  Strom 
Civil  Engineer 
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Seismic  Performance  Evaluation  of 
St.  Stephens  Hydropower  Plant 


1  Introduction 


Scope  of  Work 

At  the  request  of  the  U.S.  Army  Engineer  District,  Charleston,  the  U.S.  Army 
Engineer  Research  and  Development  Center,  Waterways  Experiment  Station  conducted 
the  seismic  performance  evaluation  of  St.  Stephen  Powerhouse  located  in  the  Cooper 
River  Rediversion  Project  in  South  Carolina  (SC).  The  structure  is  located  about  60  km 
north  of  Charleston,  SC,  and  consists  of  a  powerhouse  reinforced  concrete  structure 
founded  on  rock  and  flanked  by  rolled-fill  earth  embankments.  This  report  summarizes 
the  seismic  evaluation  of  the  erection  bay  superstructure. 

A  performance-based  approach  is  used  for  the  seismic  evaluation  to  determine  the 
risk  of  superstructure  collapse  during  a  major  earthquake.  In  the  performance-based 
approach,  superstructure  displacement  ductility  capacities  are  compared  to  earthquake 
displacement  ductility  demands.  Collapse  prevention  performance  requires  that  collapse 
of  the  structure  be  prevented  regardless  of  the  level  of  damage.  Resistance  can  decrease 
with  increasing  displacements  provided  the  structure  will  not  collapse  when  subjected  to 
extreme  earthquake  events.  The  powerhouse  superstructure  is  assumed  to  perform 
satisfactorily  (collapse  prevented)  if  member  flexural  displacement  ductility  capacities 
are  greater  than  displacement  ductility  demands,  and  if  shear  capacities  are  greater  than 
shear  demands. 

The  erection  bay  is  more  vulnerable  to  earthquake  ground  motions  than  are  the 
generator  bays.  The  erection  bay  is  assumed  to  be  critical  because  of  the  additional 
hazard  associated  with  the  possibility  of  derailment  of  the  overhead  crane.  The  overhead 
crane  will  likely  be  parked  in  the  erection  bay  at  the  time  of  a  major  earthquake  event. 

Displacements  and  displacement  ductility  demands  are  determined  by  elastic 
response  spectrum  analysis.  FEMA  273,  “NEHRP  Guidelines  for  the  Seismic 
Rehabilitation  of  Buildings”  considers  this  to  be  a  Linear  Dynamic  Procedure  (LDP)  type 
analysis.  Earthquakes  considered  in  the  collapse  prevention  analysis  are  the  maximum 
considered  earthquake,  a  probabilistic  event  with  a  return  period  of  2475-years  and  the 
Maximum  Credible  Earthquake  (MCE),  a  deterministic  event.  The  “mean  +  one  standard 
deviation”  event  was  used  for  the  MCE  analysis. 

Elastic  response  spectrum  analysis  provide  both  force  and  displacement  demands.  A 
SAP2000  finite  element  method  analysis  provided  the  elastic  seismic  demands  for  the 
two  design  response  spectra. 
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General  Description  of  the  Erection  Bay 

A  general  layout  of  the  Powerhouse  is  presented  in  Figure  1-1  and  1-2.  As  seen  in 
these  figures,  St.  Stephen  Powerhouse  consisted  of  three  generator  bays  and  one  erection 
bay.  Figures  1  shows  the  in-plan  layout  of  the  Power  plant  at  EL.  2.0.  In  this  figure, 
looking  upstream,  the  66  ft  wide  erection  bay  is  seen  on  the  left  of  the  generation  bays. 
Figure  1-2  shows  a  longitudinal  section  of  the  Powerhouse,  taken  at  a  section  along  the 
powerhouse  centerline.  This  figure  shows  the  elevations  of  different  floor  levels  and  the 
general  layout  of  the  erection  bay.  A  seismic  performance  evaluation  of  the  erection  bay 
is  conducted  in  this  investigation. 

Figure  1-3  is  a  transverse  section  through  the  middle  of  the  erection  bay.  As  seen  in 
this  figure,  the  main  erection  bay  has  several  floor  levels  over  a  massive  concrete  base 
that  extends  from  EL  1 8.0  to  EL.  -41 .  Precast  prestressed  tee-beams  at  El  93.92 
comprise  the  roof  system.  Field  observation  reported  that  the  actual  tee-beam  dimensions 
agree  with  the  PCI  standard  Type  “B”  load  table  properties.  Therefore,  the  standard  PCI 
“8ST32”  tee-beam  is  used  for  the  seismic  analysis.  An  8-ft  wide  “8ST32”  has  a  depth  of 
32  inches,  a  unit  weight  of  560  pounds  per  linear  feet,  a  section  modulus  of  4,667  in.^  for 
the  top  and  1,619  in.^  for  the  bottom,  and  a  cross-sectional  area  of  506  in.^. 

Figures  1-4  and  1-5  show  the  inside  views  of  the  Powerhouse  from  the  erection-floor 
at  El  57.0.  A  135/20-ton  overhead  crane,  manufactured  by  Broadline  Corporation, 
Richmond,  CA,  is  supported  by  corbels  located  on  the  face  of  the  upstream  and 
downstream  walls  (El.  81.25).  Internal  documents  indicate  that  the  overhead  crane  used 
at  the  site  weighs  1 52,200  pounds.  The  corbel  seat-width  must  be  large  enough  to 
prevent  the  crane  from  dropping  to  the  generator  floor  (wide  enough  to  prevent 
unseating)  during  a  major  earthquake  event.  The  crane  travels  on  rails  mounted  14-inches 
from  the  corbel  face.  Displacement  demands  on  the  wall  at  the  corbel  location  therefore 
must  be  less  than  14-inches  to  prevent  unseating. 

Photographs  taken  of  the  upstream  and  downstream  faces  of  the  powerhouse  are 
shown  in  Figures  1-6  and  1-7  respectively. 
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Figure  1-1  Cooper  River  Rediversion  project  Power  plant  -  general  arrangement  (Plan-  El  2.0) 
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Figure  1-3  Cooper  River  Rediversion  project  Power  Plant  -  Transverse  section  through  the  middle  of  the  Erection  Bay 


Figure  1-6.  Downstream  view  of  the  St.  Stephen  Powerhouse  looking  from  the  access 
road. 
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Figure  1-7.  Upstream  view  of  the  St.  Stephen  Powerhouse  looking  from  the  fish-lift  side. 
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Figure  1-8.  Upstream  view  of  the  intake-deck  looking  from  the  fish-lift  side. 
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structure  Idealization  for  Seismic  Evaluation 


In  many  cases  the  seismic  analysis  can  be  based  on  a  “superstructure  only”  analytical 
model.  The  “superstructure  only’  model  assumes  the  superstructure  acts  as  a  flexible 
appendage  attached  to  a  rigid  substructure.  When  subject  to  earthquake  ground  motions 
the  substructure,  since  it  is  considered  to  be  rigid,  moves  in  unison  with  the  ground  and 
the  accelerations  at  the  top  of  the  substructure  are  equal  to  the  accelerations  at  the 
substructure-foundation  interface.  The  flexible  superstructure  appendage  is  therefore 
assumed  to  respond  as  if  it  were  resting  directly  on  the  ground.  Although  not  completely 
rigid,  the  fundamental  mode  of  vibration  for  the  substructure  is  usually  less  than 
0.05  seconds. 

The  rigid  substructure  /  flexible  superstructure  assumption  is  expected  to  provide 
reasonable  estimates  of  the  earthquake  demands  for  the  superstructure  members  provided 
the  periods  of  vibration  for  the  modes  contributing  to  superstructure  response  are 
significantly  longer  than  the  fundamental  period  of  the  substructure.  Magnification 
effects  can  be  significant  when  any  of  the  superstructure  contributing  modes  have  periods 
of  vibration  that  are  less  than  twice  that  of  the  substructures  fundamental  mode.  With 
respect  to  the  St.  Stephen  powerhouse,  the  superstructure  is  considered  to  be  that  part  of 
the  upstream  wall  and  roof  system  above  the  intake  deck  (above  elevation  86.0)  on  the 
upstream  side,  and  that  part  of  the  downstream  wall,  attached  structures,  and  roof  system 
above  the  tailrace  deck  (above  elevation  57.0)  on  the  downstream  side.  The  short  stiff 
upstream  wall  above  elevation  86.0  makes  the  superstructure  vulnerable  to  force  and 
displacement  demand  magnification  effects.  Therefore,  a  finite  element  analytical  model 
including  the  substructure  was  necessary  to  obtain  reasonable  estimates  of  force  and 
displacement  demands  on  vulnerable  superstructure  components,  i.e.  walls. 

Usually  the  roof  support  bearings  permit  free  rotation  at  each  end  of  the  roof  system, 
but  free  translation  is  limited  to  one  end  only  and  the  extent  of  translation  limited  by 
keepers  or  slotted  plates  to  that  required  for  temperature  expansion  and  contraction.  The 
assumptions  made  with  respect  to  the  boundary  conditions  at  each  end  of  the  roof  span 
are  critical  to  the  analysis  for  walls.  If  the  earthquake  displacement  demands  on  the  wall 
exceed  the  free  translation  capacity  of  the  bearings  the  stiffer  wall  will  begin  to  carry 
some  of  the  inertial  force  due  to  the  more  flexible  wall  mass.  In  such  circumstances  an 
analytical  model  of  the  entire  roof-wall  superstructure  system  is  required.  Both  “fixed 
against  translation”  and  “free  to  translate”  boundary  conditions  must  be  investigated  to 
obtain  maximum  earthquake  demands  on  each  powerhouse  superstructure  wall. 

Field  investigations  indicate  the  tee  beams  for  the  high  roof  are  fixed  at  the  upstream 
wall  (shim  plates  welded  to  embedded  wall  plate  using  3/16  inch  fillet  welds),  and  free  at 
the  downstream  wall  (tee  beam  bearing  plate  resting  on  neoprene  pad,  bolt  through  plate 
with  slotted  holes  in  plate).  Figures  1-9  and  1-10  show  the  upstream  and  downstream 
support  conditions  and  the  available  embedded-plate  widths  supporting  the  Tee-beam. 
Because  the  bearings  are  not  ideally  fixed  or  free,  and  because  earthquake  displacement 
demands  are  expected  to  exceed  the  free  translation  capacity  of  the  expansion  end 
bearings,  two  separate  conditions  are  explored  to  determine  the  potential  influence  on 
force  and  displacement  demands.  The  first  condition,  identified  as  “the  ideal  condition”, 
assumes  the  bearings  will  perform  as  intended  with  the  bearings  at  the  upstream  wall 
fixed  against  translation  and  the  bearings  on  the  downstream  wall  free  to  translate  (see 
Figure  1-10).  For  the  ideal  condition  the  entire  roof  mass  will  be  attracted  by  the  fixed 
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support  on  the  short  upstream  wall.  This  can  be  modeled  by  eliminating  translation 
restraint  for  the  tee  beams  at  the  free  bearing  supports  (on  top  of  downstream  wall).  The 
second  condition  identified  as  the  “non-ideal”  condition  assumes  that  the  free  bearing  has 
reached  its  translation  limit  and  both  walls  will  attract  inertial  forces  in  proportion  to  their 
lateral  stiffness.  Modeling  for  this  condition  can  be  accomplished  by  preventing  joint 
restraints  at  each  end  of  the  tee  beams.  The  condition  providing  the  greatest  force  and 
displacement  demands  (“ideal  condition”  for  upstream  wall,  “non-ideal  condition”  for 
downstream  wall)  are  used  for  final  assessments. 
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Figure  1-9.  Upstream  Tee-beam  connection  -  shim  plates  welded  to  embedded  wall  (an 
ideal  hinge  condition) 
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Figure  1-10.  Downstream  Tee-beam  connection  -  Tee  beam  bearing  plate  rests 
on  neoprene  pad,  bolt  through  plate  with  slotted  holes  (an  ideal  roller  condition). 

Note:  The  slot  is  assumed  to  permit  less  than  1-inch  of  free  travel  in  each 

direction.  The  free  travel  length  is  estimated  on  the  assumption  that  the 
washer  shown  above  completely  conceals  the  slot. 


The  capacity  of  the  fixed  bearings  to  resist  lateral  load  will  be  limited  by  the  fillet  weld 
strength.  The  shear  resistance  of  a  fillet  weld  is: 

Rn  =  <|)Fw  Aw  Where: 

The  resistance  factor,  <|)  =  0.75  (per  Table  J2.5  AISC  Load  and  Resistance  Factor  Design) 
Fw  =  0.60  Fexx  =  0.60  (70)  =  42  ksi  (per  Table  J2.5  AISC) 

Aw  =  effective  area  of  weld  =  0.707  (.1875)  =  0.13  inVinch 
R„  =  0.75  (42)  0. 13  =  4  kips  /  inch 

Assume  the  weld  length  per  bearing  =  2  sides  x  6.5  inches  -  4  x  0.5  inches  =11.0  inches 

Based  on  the  above  assumptions  the  fixed  bearings  would  have  a  capacity  to  resist  1 1  (2) 
4  =  88  kips  /  per  tee  beam  or  1 1  kips  per  foot.  It  should  be  noted  that  there  are  two  stems 
per  beam  and  therefore  two  connections  per  beam. 
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System  mass 


Earthquake  force  and  displacement  demands  are  a  function  of  system  mass. 

Therefore  the  mass  of  all  structural  and  nonstructural  components  are  needed  for  the 
development  of  a  reliable  analytical  model. 

Mass  attributed  to  the  non-structural  roofing  needs  to  be  included.  For  the 
powerhouse,  roofing  consists  of  an  elastomeric  membrane  with  2-inches  of  gravel  ballast. 
A  roofing  weight  of  40  psf  for  the  superimposed  dead  (10  psf)  and  live  loads  (30  psf) 
was  used  for  the  design.  Roof  live  load  was  not  considered  in  the  seismic  evaluation 
since  roof  live  load  would  not  likely  be  present  during  an  earthquake  event. 

The  mass  of  the  structural  elements  must  also  be  estimated  accurately  to  achieve 
reasonable  estimates  of  force  and  displacement  demand.  Those  wall  sections  that  are 
continuous  from  roof  to  supporting  floor  (without  door  openings)  will  be  required  to 
carry  the  inertial  forces  generated  by  the  mass  above  door  openings.  Assuming  a  one- 
foot  section  of  wall  is  analyzed,  the  effect  of  this  additional  mass  can  be  estimated  by 
multiplying  the  mass  of  the  one-foot  section  by  the  ratio  of  total  wall  length  (including 
openings)  to  the  total  length  of  continuous  wall  sections.  For  the  St.  Stephen  Powerhouse 
a  10%  increase  in  mass  for  the  continuous  walls  seems  reasonable  since  the  extent  of 
door  openings  is  small. 

As  stated  previously  the  mass  of  the  overhead  crane  is  also  important  to  the  analysis. 
It  is  assumed  for  the  analysis  that  the  overhead  crane  is  positioned  in  erection  bay.  The 
total  weight  of  the  135/20-ton  crane  is  152,200  pounds.  The  crane  weight  is  equally 
distributed  to  the  upstream  and  downstream  walls,  with  50%  assigned  to  the  upstream 
wall  and  50%  assigned  to  the  downstream  wall.  For  the  erection  bay  this  will  result  in  a 
weight  per  foot  of  wall  (for  both  the  upstream  and  downstream  walls)  equal  to: 

A  uniform  lengthwise  distribution  along  the  entire  erection  bay  is  acceptable  since 
as  yielding  takes  place  at  one  wall  location  the  inertial  forces  will  quickly  distribute  to 
other  non-yielding  wall  sections.  This  process  will  occur  well  before  collapse.  The  crane 
mass  is  included  in  both  the  “ideal”  and  “non-ideal”  bearing  boundary  condition  models 
described  above.  The  crane  mass  will  increase  force  and  displacement  demands  and 
increase  periods  of  vibration. 

Material  properties 

From  the  drawings  it  was  determined  that  the  design  compressive  strength  for  the 
concrete  was  3,000-psi  for  the  superstructure  and  2,500-psi  the  substructure  mass 
concrete.  Grade  40  steel  (yield  strength  equal  to  40,000  psi)  was  used  for  the 
reinforcement  of  structural  concrete  as  well  to  control  cracking  due  to  temperature  and 
shrinkage. 

The  actual  concrete  compressive  strength  used  for  seismic  evaluation  purposes  will 
be  higher  than  the  design  compressive  strength.  Practice  is  to  use  an  actual  compressive 
strength  equal  to  150%  of  the  design  compressive  strength.  The  50%  increase  accounts 
for  long-term  strength  gain  (above  the  28-day  design  strength),  and  accounts  for  the 
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probable  strength  being  greater  than  design  strength.  Concrete  shear  strength,  cracking 
moment  capacity,  and  modulus  of  elasticity  are  calculated  based  on  the  actual 
compressive  strength. 

The  modulus  of  elasticity  for  concrete  for  both  structural  and  mass  concrete  is 
approximated  by  Equation  1 . 


Ec  = 

Where: 


(1-1) 


w  =  Weight  of  concrete  in  pcf  (can  be  assumed  as  equal  to  150  pcf). 

=  Actual  strength  of  concrete 

(1.5  times  the  specified  28-day  design  compressive  strength). 

Based  on  Equation  1 ,  the  modulus  of  elasticity  for  the  concrete  based  on  actual 
strength  is  estimated  to  be  555,408  ksf. 

Yielding  is  expected  to  occur  in  the  superstructure  walls.  The  effective  moment  of 
inertia,  Ie,  of  reinforced  concrete  members  at  near  yield  conditions  can  be  significantly 
less  than  that  represented  by  the  gross  section  moment  of  inertia,  Iq.  For  powerhouse 
superstructure  walls,  the  effective  moment  of  inertia  should  be  calculated  and  used  in  the 
response  spectrum  analyses  to  assure  the  response  of  the  powerhouse  superstructure  to 
earthquake  ground  motions  is  reasonable.  The  effective  moment  of  inertia  is  an  average 
value  for  the  entire  member  and  considers  the  distribution  of  cracking  along  the  member 
length.  The  effective  moment  of  inertia  of  reinforced  concrete  structures  can  be 
estimated  based  on  the  relationship  between  the  cracking  moment  (i.e.,  the  moment 
required  to  initiate  cracking  while  ignoring  the  reinforcing  steel)  and  the  nominal  moment 
capacity  of  the  reinforced  concrete  wall  section.  The  nominal  moment  capacity  is  the 
moment  capacity  of  the  cracked  section  calculated  in  accordance  with  provisions  and 
assumptions  of  the  ACI 3 1 8  strength  design  method  before  the  application  of  any 
strength  reduction  factors.  In  the  ACI  strength  design  method  it  is  assumed  that  strains 
are  linear  across  the  section.  The  nominal  moment  capacity  for  pure  bending  assumes  that 
the  ultimate  capacity  of  the  concrete  in  compression  (compressive  stress  block)  is  equal 
in  magnitude  and  opposite  in  direction  to  the  ultimate  (yield)  capacity  of  the  reinforcing 
steel.  The  concrete  ultimate  compressive  force  and  steel  ultimate  tensile  force  create  a 
couple  that  provides  the  nominal  moment  capacity  for  the  section.  The  nominal  moments 
and  cracking  moments  used  to  estimate  effective  moment  of  inertia  are  for  those  regions 
where  moments  are  at  their  maximums.  For  powerhouse  walls  this  can  be  at  the  base  of 
cantilever  wall  members  or  in  upper  wall  sections  where  an  abrupt  change  in  wall 
thickness  occurs.  Once  the  cracking  moment  (Mcr)  and  the  nominal  moment  capacity 
(Mn)  have  been  determined  the  ratio  of  the  effective  stiffness  (Ie)  to  the  gross  stiffness 
(Iq)  can  be  estimated  as  follows: 


0.8 -0.9 


Mjl 

^CR 


Must  be  >  0.35  and  <  0.80 


(1-2) 
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Where: 


Ig  =  Gross  section  moment  of  inertia. 

Ie  =  Effective  moment  of  inertia. 

Mn  =  Nominal  moment  capacity 
Mcr  =  ^Cracking  moment  capacity. 

♦Moment  required  to  crack  a  structural  concrete  member  assuming  no  contribution  from  the  reinforcement, 
and  assuming  a  modulus  of  rupture  of  7.5'rJ/^ 

Sectional  properties  for  the  erection  bay  wall  elements 

The  following  material  properties  have  been  determined  for  the  erection  bay  walls 
assuming  that  the  walls  will  have  earthquake  moment  demands  in  excess  of  moment 
capacity,  i.e.  yielding  occurs: 

Cracking  moment  capacity 
Nominal  moment  capacity 

Ratio  of  effective  moment  of  inertia  to  gross  moment  of  inertia  (Ie/Iq) 


These  sectional  properties  are  summarized  in  Table  1-1  based  on  actual  compressive 
strength  conditions.  For  the  initial  SAP  2000  analysis  all  the  members  not  listed  in  Tables 
1-1  or  1-2  will  be  assigned  an  effective  stiffness  equal  to  0.8  times  the  gross  stiffness  (EIe 
=  0.8  EIg).  The  use  of  a  maximum  effective  stiffness  equal  to  80-percent  of  the  gross 
section  stiffness  assumes  that  some  cracking  and  slight  degradation  in  stiffness  will  occur 
under  cyclic  earthquake  loading.  The  20-percent  minimum  reduction  in  stiffness  is  per 
recommendations  by  M.  J.  N.  Priestley.  Since  the  period  of  vibration  is  equal  to  the 
square  root  of  the  mass  divided  by  the  stiffness  it  can  be  concluded  that  a  20-percent 
reduction  in  stiffness  will  provide  earthquake  demands  that  vary  only  slightly  from  those 
obtained  from  analyses  based  on  gross  section  properties.  The  precast  /  prestressed 
concrete  tee  beams  should  not  crack  during  an  earthquake.  Therefore  for  those  members 
the  gross  section  properties  described  in  Table  2  can  be  used  in  the  analysis.  The 
nominal  moment  capacities  are  based  on  the  ultimate  capacity  of  the  concrete  at  a 
maximum  strain  of  0.002,  and  on  the  yield  strength  of  the  steel  assuming  an  ultimate 
strain  capacity  of  5-percent.  By  limiting  the  compressive  strain  in  the  concrete  to  0.002 
the  opportunity  for  compressive  micro  cracking,  which  can  lead  to  bond  deterioration  and 
eventually  to  rebar  splice  failure  in  plastic  hinge  regions,  is  prevented. 
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Table  1-1 

Sectional  properties  for  wall  elements  -  fca  =  4500  psi 
(Actual  compressive  strength) 


Wall 

Thickness 

(Feet) 

Nominal 

Moment 

Capacity  -  Mn 
(Foot-Kips) 

Cracking 

Moment 

Capacity  -  Mcr 
(Foot-Kips) 

Ie/Ig 

1.0 

17.23 

12.83 

0.49 

1.5 

32.78 

30.04 

0.72 

3.5 

161.60 

169.60 

0.80 

4.5 

225.20 

278.30 

0.80 

The  properties  for  the  tee  beams  as  described  in  Table  2  are  used  for  the  analysis. 
The  mass  of  the  tee  beams  can  be  uniformly  distributed  along  the  tee  beam  length,  or  can 
be  modeled  as  a  series  of  lumped  masses. 


Table  1-2 

Tee  Beam  Properties 


Beam  Type 

Area  (A) 

(in') 

Moment 
of  Inertia  (I) 

(in") 

Modulus  of 
Elasticity  (E) 

(ksi) 

wt. 

(Kips/  foot) 

32”  single  tee 

(high  roof) 

538. 

49,329. 

4,287. 

0.560 

16”  double  tee 

(low  roof) 

325. 

6,634. 

4,287. 

0.339 

*  Tee  beam  properties  are  for  an  8.0’  wide  section. 
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2  Performance-Based  Evaluation 


General 

A  performance-based  approach  is  used  to  determine  if  the  St.  Stephen  powerhouse 
erection  bay  will  meet  collapse  prevention  performance  objectives  when  subject  to 
maximum  considered  earthquake  (2475-year  event),  and  Maximum  Credible  Earthquake 
(MCE)  ground  motion  demands.  The  procedure  requires  computation  of  displacement 
and  displacement  ductility  demands  for  vulnerable  powerhouse  members,  i.e., 
superstructure  walls. 

Performance  objectives  are  met  if: 

•  Brittle  modes  of  failure  are  suppressed  to  allow  flexural  yielding  to  take  place  in 
the  superstructure  walls. 

•  Member  displacement  ductility  capacity  exceeds  displacement  ductility 
demands. 


Detailed  discussions  on  the  performance-based  seismic  analysis  and  its  applications 
in  seismic  assessment  for  building  structures  are  available  in  the  literature  (Medhekar  and 
Kennedy  2000,  Chandler  and  Mendis  2000,  and  Moehle  1992). 


Procedures  Used  to  Assess  Structure  Displacement 
Demand 

A  Linear  Dynamic  Procedure  (LDP)  analysis  for  the  St.  Stephen  powerhouse  will  be 
used  to  assess  seismic  performance.  The  LDP  uses  a  response  spectrum  analysis  to 
determine  force  and  displacement  demands  on  the  structure.  The  LDP  is  described  in 
detail  in  FEMA  356,  “Prestandard  and  Commentary  for  the  Seismie  Rehabilitation  of 
Buildings”.  The  desired  performance  level  for  this  evaluation  is  “collapse  prevention”. 
Collapse  prevention  permits  yielding  to  occur,  provided  such  yielding  does  not  lead  to 
collapse  of  the  structure.  Displacements  less  than  those  causing  collapse  are  required  as 
illustrated  in  Figure  2-1. 
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Increasing  earthquake  demand 


> 


Figure  2-1  Performance  and  Structural  Deformation  Demand  for  Ductile  Structures 


Earthquake  Demands  on  Inelastic  Systems 

A  Linear  Dynamic  Procedure  (LDP),  i.e.,  response  spectrum  analysis  is  a  linear 
elastic  analysis  that  provides  force  and  displacement  demands  assuming  the  structure 
remains  elastic  during  the  design  earthquake.  The  inelastic  response  of  a  structure 
subjected  to  earthquake  ground  motion  however  is  different  than  the  elastic  response.  The 
difference  occurs  because  the  vibrational  characteristics  of  the  structure  change  as  the 
structure  yields.  The  predominant  change  is  a  shift  in  the  fundamental  period  of 
vibration.  In  most  cases,  a  reduction  in  earthquake  demand  occurs  as  the  period  of  the 
structure  lengthens.  In  Figure  2-2,  a  capacity  spectrum  is  used  to  illustrate  the  inertial 
force  reduction  (or  spectral  acceleration  reduction)  that  occurs  when  a  structure  yields.  In 
Figure  2-2,  earthquake  demands  for  an  elastic  system,  as  represented  by  a  response 
spectrum,  are  reconciled  with  the  elasto-plastic  load/displacement  characteristics  for  a 
ductile  structure.  Point  "A"  represents  the  earthquake  demand  assuming  the  structure 
remains  elastic  with  the  line  "0-A"  the  linear  elastic  response.  Point  "B"  represents  the 
earthquake  demands  for  elasto-plastic  behavior  with  the  line  "O-B-C"  representing  the 
load/displacement  characteristics  of  the  elasto-plastic  system. 

As  can  be  seen,  the  earthquake  demand  on  the  elasto-plastic  system  is  substantially 
less  than  that  of  the  elastic  system.  However,  when  the  fundamental  period  of  the 
structure  falls  in  the  ascending  portion  of  the  response  spectrum,  a  shift  in  period  may 
actually  increase  earthquake  demands.  This  can  occur  for  stiff  structures  founded  on  soft 
soils.  It  should  also  be  recognized  that  damping  increases  significantly  as  structures  yield. 
This  also  reduces  earthquake  demand.  Simple  techniques  have  been  developed  for 
estimating  earthquake  demands  on  inelastic  systems.  For  structures  with  fundamental 
periods  of  vibration  located  in  the  constant  acceleration  range  of  the  response  spectrum, 
the  inelastic  response  capability  will  be  a  function  of  structure  displacement  ductility  (|li). 
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and  can  be  estimated  using  the  equal  energy  principle.  The  fundamental  period  of 
vibration  of  the  structure  is  designated  as  “T”  in  this  report.  For  structures  with 
fundamental  periods  of  vibration  in  the  constant  velocity  range  of  the  response  spectrum, 
the  inelastic  response  capability  will  also  be  a  function  of  structure  displacement  ductility 
(p),  and  can  be  estimated  using  the  equal  displacement  principle. 


Figure  2-2  Earthquake  Demands  on  Inelastic  Structures 


Structures  with  fundamental  periods  of  vibration  less  than  the  characteristic  ground 
motion  period  will  exhibit  an  equal  energy  response.  The  characteristic  ground  motion 
period  is  designated  as  “To”  in  this  report  and  represents  the  intersection  of  the  response 
spectrum  constant  acceleration  range  and  constant  velocity  range.  The  characteristic 
ground  motion  period  generally  varies  between  0.2  seconds  and  0.7  seconds  depending 
on  site  conditions,  with  firm  sites  having  shorter  characteristic  periods  than  soft  sites. 
The  characteristic  ground  motion  period  is  illustrated  in  Figure  2-3.  This  figure  is  a 
tripartite  plot  of  a  given  site  specific  response  spectrum  which  has  been  smoothed  to 
provide  a  straight  line  (on  log  /log  scale)  representation  of  earthquake  demands  in  the 
constant  acceleration,  constant  velocity,  and  constant  displacement  ranges. 
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Figure  2-3  Characteristic  Ground  Motion 

For  an  equal  energy  response,  the  structure  must  have  sufficient  displacement 
ductility  to  provide  the  reserve  inelastic  energy  capacity  needed  to  resist  earthquake 
ground  motion  demands.  The  equal  energy  response  concept  is  presented  in  Figure  2-4. 
For  a  given  displacement  ductility  (jxa)  the  inelastic  (yield)  capacity  (Fy)  must  be 
sufficient  to  produce  an  equal  energy  response.  Equating  the  energy  for  a  linear  elastic 
response  to  that  for  an  inelastic  response  (hatched  area  under  the  nonlinear  portion  of  the 
load  displacement  curve  equal  to  the  hatched  area  under  the  linear  elastic  curve),  it  can  be 
determined  that  the  yield  capacity  of  the  structure  must  be  equal  to  or  greater  than  the 

capacity  required  of  the  structure  if  it  were  to  remain  elastic  (Fe)  divided  by  - 1  , 

or: 
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Figure  2-4.  Equal  Energy  Response 

It  can  also  be  seen  that  the  displacement  capaeity  of  the  inelastic  structure  must  be  larger 
than  the  displacement  of  the  elastic  structure  if  the  structure  is  to  survive. 

Structures  with  fundamental  periods  of  vibration  greater  than  the  characteristic  ground 
motion  period  will  exhibit  an  equal  displacement  response.  An  equal  displacement 
response  means  that  to  perform  as  intended,  the  displaeement  ductility  capacity  must  be 
sufficient  to  provide  a  structure  displacement  capacity  equal  to,  or  greater  than,  the  peak 
displacement  the  structure  will  experience  during  the  design  earthquake.  The  equal 
displacement  response  concept  is  presented  in  Figure  2-5.  From  Figure  2-5  it  ean  be 
determined  that  the  yield  capacity  of  the  structure  must  be  equal  to  or  greater  than  the 
capacity  required  of  the  structure  if  it  were  to  remain  elastic  (Fe)  divided  by  ps,  or: 
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Figure  2-5.  Equal  Displacement  Response 


It  can  also  be  seen  that  the  displacement  of  the  inelastic  structure  equals  that  of  the  elastic 
structure. 


A  general  relationship  has  been  developed  (Paulay  and  Priestley,  1992)  for  relating 
required  yield  strength  to  elastic  demand.  This  relationship  provides  a  smooth  transition 


from  an  equal  acceleration  response 


( 


Fy 


=  1  regardless  of  ps)  at  T  =  0,  through  the  equal 


energy  approximation  (— ^  = 


- 1 )  at  about  T  =  0.75  To,  to  the  equal  displacement 


approximation  (— ^  =  //^ )  for  T  >  1 .5  Tq.  The  relationship  for  the  smooth  transition  is: 


Fy  i.sr, 


(2-3) 
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A  form  of  the  above  equation  was  used  to  determine  the  flexural  displacement  ductility 
demand  (ps)  on  the  St.  Stephen  erection  bay  walls.  The  equation  is  designated  in  the 
calculations  as  the  Priestley  Equation.  Since  flexural  demands  are  the  quantities  of 
interest  when  determining  flexural  displacement  ductility  demand,  the  above  equation  is 
expressed  in  terms  the  moment  demand  obtained  from  a  linear  elastic  response  spectrum 
(modal)  analysis  (Mde)  and  the  yield  or  nominal  moment  capacity  of  the  particular 
member  under  investigation  (Mn).  Therefore  for  a  flexural  response,  the  relationship  for 
the  smooth  transition  is  in  terms  of  the  displacement  ductility  demand  (ps)  is: 


Us  =1.5 


M 


DE 


M. 


T 


(2-4) 


The  above  equation  was  used  to  estimate  the  displacement  ductility  demands  on  critical 
structural  elements  of  the  St.  Stephen  Powerhouse  erection  bay. 


Ductility  Capacity  of  Reinforced  Concrete  Structures 

Performance  evaluations  require  displacement-based  analyses  to  determine  if 
member  flexural  displacement  ductility  capacities  are  greater  than  flexural  displacement 
ductility  demands.  Flexural  displacement  ductility  capacity  is  related  to  curvature 
capacity,  plastic  hinge  length  (length  of  zone  where  yielding  occurs),  and  member  length. 

The  curvature  capacity  will  depend  on  the  maximum  amount  of  strain  that  can  be 
placed  on  the  concrete  and  reinforcing  steel.  To  prevent  fracturing  of  the  reinforcing 
steel,  the  steel  strain  should  be  limited  to  5-percent.  Static  load  testing  produces  ultimate 
strain  capacities  of  18-percent  or  more,  but  due  to  low  cycle  fatigue  effects  the  maximum 
permitted  for  earthquake  load  conditions  is  5-percent.  The  maximum  ultimate  strain 
capacity  of  concrete  under  dynamic  loadings  is  approximately  0.4-percent.  However, 
when  there  is  inadequate  confinement  steel,  the  maximum  ultimate  strain  capacity  should 
be  limited  to  0.2-percent  to  prevent  compressive  micro  cracking  of  the  concrete  that  can 
lead  to  bond  strength  deterioration  and  to  subsequent  rebar  splice  failure.  Using  the  above 
strain  limit  information,  the  ultimate  curvature  capacity  of  a  reinforced  concrete  section 
can  be  determined  (see  Figure  2-6). 
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Figure  2-6.  Ultimate  Curvature  Capacity 


Ultimate  rotation  capacity  is  equal  to  the  ultimate  curvature  capacity  times  the  plastic 
hinge  length.  When  the  nominal  moment  capacity  (Mn)  is  less  than  1.2  times  the  cracking 
moment  (Mcr),  the  plastic  hinge  length  to  be  used  in  calculating  rotational  capacity  is: 

Ip  =QM7)5fy{d^)  (Mpa  units)  (2-5) 

|/^  =  0.30 fy  (df,  )|  (Ksi  units),  where: 

fy  =  yield  strength  of  the  reinforcing  steel 
db  =  diameter  of  reinforcing  steel 

When  the  nominal  moment  (Mn)  is  greater  than  twice  the  cracking  moment  (Mcr),  the 
plastic  hinge  length  to  be  used  in  calculating  rotational  capacity  is: 

/^  =  0.08  Z -I- 0.021 8/y(r/*)  mm  (Mpa  units)  (2-6) 

Ip  =  0.08 Z  -I-  0.1 5 fy  (d/, )  Inches  (Ksi  units) 

Where:  L  =  length  of  the  member. 
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For  nominal  moment  strengths  between  1.2  Mcr  and  2.0  Mcr  the  plastic  hinge  length  can 
be  determined  by  linear  interpolation  between  the  results  obtained  from  the  two  plastic 
hinge  length  equations.  The  ultimate  rotational  capacity  of  the  member  can  be  estimated 
as  follows: 


(2-7) 

Where: 

0u  =  Ultimate  rotational  capacity 
(|)u  =  Ultimate  curvature  capacity 
Ip  =  Plastic  hinge  length 

For  cantilevered  structures  such  as  powerhouse  walls  it  is  advantageous  to  determine 
ultimate  displacement  ductility  capacities  and  use  those  to  compare  with  the  displacement 
ductility  demands  obtained  from  a  seismic  analysis. 

The  displacement  ductility  capacity  of  powerhouse  superstructure  walls  and  other 
members  that  behave  as  cantilever  structures  can  be  determined  using  various  cantilever 
beam  models.  For  structures  where  the  mass  can  be  assumed  to  occur  at  the  top  of  the 
structure  (concentrated  mass  model)  Figure  2-7  applies.  This  model  was  used  to  evaluate 
the  displacement  ductility  capacity  of  the  St.  Stephen  powerhouse  walls. 


8u 


Figure  2-7.  Concentrated  Mass  Model 


For  the  concentrated  mass  model  (Figure  2-7)  the  ultimate  displacement  capacity  can  be 
estimated  by: 
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(2-8) 


Or  the  ultimate  displacement  ductility  capacity  can  be  estimated  by: 


//  —  1  -I"  3 


^-1 

(  1 
1-0.5-^ 

) 

1 

1  n 

(2-9) 


Where: 

(|)y  =  Curvature  at  first  yield  of  the  reinforcing  steel. 

1  =  Length  of  cantilever 

The  above  equations  were  used  to  assess  the  displacement  ductility  capacity  and 
displacement  ductility  capacity  of  the  St.  Stephen  erection  bay  walls. 


Moment  -  Curvature  (M-Phi)  Analysis 

The  maximum  curvature  capacity  is  needed  to  determine  the  maximum  displacement 
ductility  capacity  for  critical  structural  elements  of  the  St.  Stephen  Powerhouse  erection 
bay.  Maximum  curvature  capacity  is  based  on  the  strain  limits  established  above  for  the 
concrete  and  steel,  and  is  dependent  on  the  axial  load  on  the  member. 

The  maximum  displacement  ductility  demand  is  a  function  of  the  nominal  moment 
capacity  of  the  member  as  indicated  by  Equation  2-4. 

Both  the  maximum  curvature  capacity  and  nominal  moment  capacity  can  be 
determined  by  moment  curvature  analysis.  Moment  curvature  analyses  for  critical 
structural  elements  of  the  St.  Stephen  Powerhouse  erection  bay  were  performed  using  the 
Corps  M-Phi  computer  program  developed  by  Mohammad  R.  Ehsani  and  Matthew  E. 
Marine.  In  a  moment  curvature  analysis,  as  with  most  section  analyses  the  strain  is 
assumed  to  vary  linearly  across  the  cross  section  of  the  member.  Strain  values  for 
concrete  and  steel  are  selected,  and  based  on  known  relationships  between  stress  and 
strain  the  forces  in  the  concrete  and  steel  are  calculated  and  resisting  moments 
determined  for  a  particular  axial  load  condition  specified  for  the  member.  The  strain 
values  selected  determine  the  curvature  as  illustrated  in  Figure  2-6.  This  process  is 
repeated  until  the  ultimate  strain  capacity  of  either  the  steel  or  concrete  is  reached  and 
until  a  complete  moment-curvature  relationship  for  the  member  is  obtained.  The  results 
of  the  M-Phi  analyses  for  the  St.  Stephen  Powerhouse  critical  members  are  illustrated  in 
Figure  2-8.  As  shown  in  that  figure,  there  is  a  reduction  in  the  moment  values 
immediately  following  cracking  of  the  section.  This  occurs  as  the  curvature  increases. 
This  behavior  is  common  for  lightly  reinforced  sections.  As  soon  as  the  concrete  cracks 
in  tension,  the  share  of  the  tensile  stress  carried  by  the  concrete  shifts  to  the  steel 
resulting  in  a  lower  moment.  Once  cracking  occurs  and  as  curvature  increases  all  the  steel 
eventually  reaches  strains  greater  than  yield.  For  lightly  reinforced  wall  sections  the  steel 
will  reach  ultimate  strain  capacity  before  the  concrete. 
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Figure  2-8  Results  of  M-Phi  analyses  for  St.  Stephen  erection  bay  walls 


Interpreting  Analysis  Results 

It  is  expected  that  flexural  yielding  will  occur  in  various  members  of  the  powerhouse 
superstructure  during  major  earthquake  events.  However,  for  flexural  yielding  to  occur 
various  brittle  modes  of  failure  such  as  shear  must  be  suppressed.  The  displacement 
ductility  demand  and  capacity  information  provided  above  relates  to  the  flexural  yielding 
mode  of  failure.  The  location  where  yielding  occurs  (plastic  hinge  regions)  can  be  easily 
identified,  being  those  locations  where  elastic  moment  demands  exceed  nominal  moment 
capacities.  For  the  MCE  +  a  and  the  2475-year  earthquake  events  flexural  yielding 
occurred  in  the  powerhouse  erection  bay  superstructure  walls  at  the  locations  where: 
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•  The  uppermost  upstream  wall  ( 1 .5-foot  thick  wall  supporting  the  roof  tee  beams) 
connects  to  the  intake  deck  at  Elevation  86. 

•  The  uppermost  downstream  wall  ( 1 .5-foot  thick  wall  supporting  the  roof  tee  beams) 
connects  to  the  crane  rail  corbel  at  Elevation  81.25. 

•  The  lower  section  of  the  downstream  wall  (3.5-foot  thick  wall)  attaches  to  the  tailrace 
deck  at  Elevation  57. 

Failure  at  any  one  of  these  yield  regions  can  lead  to  collapse  of  the  erection  bay  roof. 
However,  to  determine  if  collapse  might  occur  during  extreme  earthquake  event,  several 
brittle  modes  of  failure,  in  addition  to  the  flexural  mode  of  failure,  must  be  also 
investigated.  The  potential  for  failure  is  assessed  using  demand  to  capacity  evaluations. 
The  type  of  demand  to  capacity  evaluations  used  to  evaluate  collapse  performance,  and 
maximum  permissible  demand  to  capacity  ratio  (OCR)  are  presented  in  Table  2-1. 

Table  2-1  Demand  to  capacity  evaluations 


Potential  Failure 
Mechanism 

Type  of  Demand  to  Capacity 
Evaluation 

Max. 

DCR 

Fracture  of  reinforcing 
steel 

Displacement  ductility  demand  to 
displacement  ductility  capacity 

1.0 

Crushing  and  spalling 
of  concrete 

Displacement  ductility  demand  to 
displacement  ductility  capacity 

1.0 

Failure  of  rebar  splices 

Displacement  ductility  demand  to 
displacement  ductility  capacity 

1.0 

Shear  failure 

Elastic  shear  demand  to  shear  capacity 

1.0 

Sliding  shear  failure 

Elastic  shear  demand  to  sliding  shear 
capacity 

1.0 

As  with  older  structures,  the  reinforcing  steel  in  the  walls  of  the  St.  Stephen 
powerhouse  is  spliced  at  the  base  of  the  walls.  This  means  the  splices  are  located  in 
potential  plastic  hinge  (yield)  regions.  Splice  failures  can  occur  as  damage  accumulates  in 
plastic  hinge  regions.  In  general  some  limited  flexural  yielding  can  be  tolerated  in 
conventionally  reinforced  structures.  But  to  assure  the  structure  will  not  collapse  due  to 
splice  failures  the  compressive  strain  in  the  concrete  confining  the  splice  should  not 
exceed  0.2-percent.  Strains  below  0.2-percent  will  prevent  micro  cracking  in  the 
concrete,  and  as  such  will  prevent  the  onset  of  bond  deterioration,  i.e.  the  splitting  of  the 
concrete  surrounding  the  steel,  that  can  lead  to  splice  failure.  Crushing  of  the  concrete 
will  not  occur  until  compressive  strains  in  the  concrete  exceed  0.4-percent.  Therefore  the 
0.2-percent  strain  limit  placed  on  the  concrete  to  prevent  splice  failures  will  also  protect 
against  spalling. 
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Shear  failures  are  brittle  sudden  failures  and  therefore  are  unacceptable.  To  prevent 
shear  and  sliding  shear  failures  the  elastic  shear  demands  from  the  response  spectrum 
analysis  should  not  exceed  the  shear  or  sliding  shear  capacity  of  the  member,  (DCR  < 

1 .0).  Shear  and  sliding  shear  capacities  are  described  below. 


Shear  Capacity 


The  capacity  of  the  concrete  in  shear  (diagonal  tension)  is  due  to  aggregate 
interlock,  to  a  lesser  extent  due  to  the  shear  resistance  available  from  the  transverse  rein¬ 
forcing  (traditional  truss  mechanism)  and  due  to  shear  strength  enhancement  as  the  result 
of  axial  load.  The  total  ultimate  shear  strength  (Fy)  can  be  taken  as: 

V^=^(Vc+Vs)  =  0.S5(Vc  +V,)  (2-10) 

A  conservative  representation  of  the  concrete  shear  strength  is: 


Where: 


k  + 


2000A„ 


(psi  units) 


(2-11) 


P  =Axial  load  on  section 

/J  =  Actual  concrete  compressive  strength  (The  actual  concrete  compressive 
strength,  which  may  be  as  high,  or  higher  than  1.5  times  the  design  compressive  strength, 
should  be  used  when  calculating  the  shear  capacity.) 

Ag  =  Gross  concrete  area.  Ae  =  0.8(Ag) 

Where  k  =  1  for  pa  =  U  and  k=  0.5  for  ps  =  2.0,  with  linear  interpolation  between  these 
values  for  ps  greater  than  1 .0  but  less  than  2.0,  where: 

Ps  =  Flexural  displacement  ductility  demand. 

The  relationship  between  concrete  shear  strength  and  flexural  displacement  ductility  is 
illustrated  in  Figure  2-9.  For  the  St.  Stephen  Powerhouse  the  lowest  shear  strength 
(residual  shear  strength)  was  used  in  the  DCR  analysis  since  the  flexural  ductility 
demands  always  exceeded  two. 
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(psi  units) 


Flexural  Displacement  Ductility  Demand  (p) 
Figure  2-9  Concrete  Shear  Strength  vs.  Ductility 


Sliding  shear  (shear  friction)  along  the  base  of  a  structure  or  structural  member 
should  also  be  investigated.  The  shear  friction  shear  capacity  (Vsf)  can  be  determined  by 
the  following  expression; 

VSF  =  ^SF  (P  +  0.25A3fy)  (2-12) 


Where: 


PsF  =  sliding  shear  coefficient  of  friction,  per  ACI 315. 

P  =  Axial  load  on  section. 

As  =  Area  of  the  longitudinal  reinforcing  steel  across  the  potential  failure  plane, 
fy  =  yield  strength  of  the  reinforcing  steel 
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3  Earthquake  Ground  Motions 


Earthquake  ground  motion  response  spectra  representing  the  probabilistically  derived 
2475-year,  1000-year,  and  475-  year  events  (exposure  periods)  were  developed  to 
represent  earthquake  demands  to  be  considered  for  the  seismic  evaluation  of  the  St. 
Stephen  powerhouse.  Deterministically  derived  Maximum  Credible  Earthquake  (MCE) 

mean,  and  mean  +  IG  (probability  of  non-exceedance  of  84%),  events  were  also 
developed  as  potential  earthquake  demands  for  the  seismic  evaluation.  These  design 
earthquake  ground  motions  were  developed  for  this  analysis  by  a  review  of  the  current 
understanding  of  the  seismic  setting  and  earthquake  sources  that  dominate  the  seismic 
hazard  for  this  site.  Earthquake  demands  are  provided  in  terms  of  acceleration  response 
spectra  for  site  class  B  conditions.  Site  class  B  represents  rock  with  a  shear  wave 
velocity  between  2500  and  5000  feet  per  second.  Acceleration  response  spectra  ordinates 
for  the  aforementioned  returned  periods  are  presented  in  Table  3-la.  Since  performance 
evaluation  was  for  collapse  prevention  only  the  2475-year  event  information  was  used  in 
the  probabilistic  response  spectrum  analyses. 

The  project’s  seismic  hazard  is  dominated  by  and  located  within  the  Charleston 
seismic  zone.  A  source  zone  is  specified  to  quantify  the  recognized  seismic  hazard 
because  the  central  and  eastern  North  American  geologic  setting  and  current  data  do  not 
allow  specifying  definitive  locations  of  causal  fault  sources.  The  Charleston  seismic 
zone  maximum  earthquake  hazard  is  predominately  based  on  the  historic  Charleston 
earthquake  of  1886  with  an  estimated  moment  magnitude  (Mw)  =  7.3  with  a  return 
period  of  650  years.  The  estimated  location  was  near  Summerville,  SC.  The  study  of 
historic  and  recorded  seismicity,  paleoseismic  studies,  and  consideration  of  the 
uncertainties  in  the  causal  seismic  sources  current  zone  boundaries  were  used  to  develop 
the  probabilistically  based  1997  National  Seismic  Hazard  maps  (Frankel  et  al,  1996).  The 
resulting  hazard  map  is  shown  in  Figure  3.1  for  PGA  with  a  return  period  of  2475  years. 

Once  the  International  Building  Code  (IBC  2000)  provisions  are  adopted  earthquake 
design  ground  motions  will  be  based  on  the  maximum  considered  earthquake.  These 
maps  are  a  derivative  of  the  NEHRP  1997  maps  for  a  2475-year  exposure,  modified 
deterministically  near  active  faults  of  the  western  US.  Because  this  site  is  not  located  near 
a  mapped  active  fault  the  maximum  considered  earthquake  maps  provide  the  same 
motions  as  the  2475  year  NEHRP  maps  for  mapped  0.2  and  1.0  second  periods. 

Maximum  considered  earthquake  demands  are  provided  in  terms  of  acceleration  response 
spectra  for  site  class  A  and  B  conditions.  Site  class  A  represents  a  hard  rock  site  with  a 
shear  wave  velocities  greater  than  5000  feet  per  second.  Acceleration  response  spectra 
ordinates  for  the  maximum  considered  earthquake  are  presented  in  Table  3- lb. 
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The  source  zone  used  to  include  the  maximum  historic  hazard  in  the  probabilistic 
maps  is  shown  in  Figure  3-2.  Within  this  broad  source  zone  there  exists  a  seismic 
“hotspot”  or  concentration  of  seismicity  near  Summerville,  SC.  This  is  the  estimated 
location  of  the  historic  1868  Charleston  earthquake.  A  recent  seismologic  evaluation  for 
this  site  (Krinitzsky  et  al,  1998)  used  this  smaller  zone  in  their  deterministic  study, 
proscribing  a  Maximum  Credible  Earthquake  (MCE)  of  7.5  Mw  occurring  within  the 
Summerville  hotspot  at  a  closest  distance  of  55  km  from  the  site.  This  deterministic 
approach  was  dictated  by  US  ACE  regulation  due  to  the  criticality  of  the  embankment 
dam  that  was  under  evaluation. 

An  updated  review  of  current  expert  opinions  for  this  area  resulted  in  a  MCE  of  7.3 
Mw  occurring  near  the  Somerville  hotspot  at  a  closest  distance  of  35  km  from  the  site. 
This  distance  was  based  on  a  zone  defined  by  the  approximate  boundary  of  the  epicentral 
region  of  the  1886  event,  Obermeier,  1999,  USCE  Contract  report.  The  source  ground 
motion  was  propagated  to  the  site  with  an  appropriate  relationship  for  CEUS,  Toro,  et  al, 
1997.  Ordinates  for  the  mean  and  mean  +1  standard  deviation  response  spectra  are 
provided  in  Table  3-1  c.  The  difference  between  these  two  recent  deterministie 
evaluations  is  the  use  of  a  0.2  Mw  units  decrease  in  the  best  estimate  of  the  1886  event 
and  the  use  of  the  closer  distance  to  the  source  boundary  of  35  km  compared  with  55  km. 
These  differences  generally  result  with  a  moderate  increase  in  the  resulting  deterministic 
hazard  that  is  in  better  agreement  with  the  probabilistie  hazard  with  a  2475-year  return 
period  that  is  proscribed  by  various  recent  building  codes  for  design  with  performance 
goals  of  eollapse  prevention. 

Aeceleration  response  spectra  ordinates  for  the  Maximum  Credible  Earthquake  are 
presented  in  Table  3-1  c.  These  deterministic  response  speetra  were  developed  to  compare 
with  2475-year  equal  hazard  spectra  developed  for  the  collapse  prevention  performanee 
evaluation.  The  comparison  indicates  the  2475-year  motions  are  a  safe  bound  and 
reasonable  when  considered  with  respect  to  the  maximum  historic  event.  A  response 
speetrum  analysis  was  also  performed  using  the  Maximum  Credible  Earthquake  “mean  + 
ct”  event. 

Site  class  B  equal  hazard  aeceleration  spectra  for  return  periods  of  475,  975,  and 
2475  years,  the  probabilistic  maximum  considered  earthquake  for  site  classes  A  and  B, 
and  the  Maximum  Credible  Earthquake  “mean”  and  “mean  +  a”  events  are  shown  in 
Figure  3-3.  These  eorrespond  to  the  information  provided  in  Tables  3-la,  3-lb,  and  3-le 
respectively. 

The  eharaeteristie  ground  motion  period,  designated  as  “To”  in  this  report,  represents 
the  intersection  of  the  response  spectrum  constant  acceleration  range  and  constant 
velocity  range.  Values  of  To  equal  to  0.20  seconds  were  selected  from  Figure  3-3  to 
represent  the  2475-year  event  and  Toro  Mean  +  a  event  (designated  as  MCE  +  a  in 
Chapter  5).  The  To  values  are  needed  for  determining  displacement  ductility  demand. 
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Figure  3-1.  NEHRP  USGS  National 
Probabilistic  Seismic  Hazard  Map  for 
PGA  with  2475  Year  Return  Period 
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Figure  3-2.  Seismicity  and  source  zones 

near  project  site 
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Figure  3-3  Earthquake  response  spectra  for  St.  Stephen  powerhouse 

Note:  For  Figure  3-3,  MCE  =  maximum  considered  earthquake.  Elsewhere  in  the  report 
“MCE”  is  used  to  represent  the  Maximum  Credible  Earthquake. 


Table  3-la  Probabilistic  acceleration  response  spectra  per  1996  NEHRP 
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Table  3-lb 


Probabilistic  acceleration  response  spectra 

Maximum  considered  earthquake  for  site  classes  A  and  B 

Return  period  =  2475-years 


Probabilistic  Acceleration  Response  Spectra 

Maximum  Considered  Earthquake  -  2475  year  event 

By  FEMA  302  (1997  NEHRP) 

Period* 

(Seconds) 

Spectral  Acceleration 

Site  Class  A  (g) 

Spectral  Acceleration 

Site  Class  B  (g) 

0.00 

0.500 

0.625 

0.06 

1.250 

1.562 

0.29 

1.250 

1.562 

0.40 

0.903 

1.129 

0.50 

0.723 

0.903 

0.70 

0.516 

0.645 

0.80 

0.452 

0.565 

1.00 

0.361 

0.452 

1.20 

0.301 

0.376 

1.50 

0.241 

0.301 

2.00 

0.181 

0.226 

Table  3-lc  Deterministic  acceleration  response  spectra 


Deterministic  Acceleration  Response  Spectra 

Toro  et  al  1997 

Period* 

Maximum  Credible 

Maximum  Credible 

(Seconds) 

Mean  +  a  (g) 

Mean  (g) 

0.00 

0.692 

0.251 

0.03 

1.540 

0.517 

0.04 

1.615 

0.559 

0.10 

1.414 

0.516 

0.20 

1.235 

0.434 

0.40 

0.938 

0.308 

1.00 

0.562 

0.183 

2.00 

0.327 

0.111 
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4  Finite  Eiement  Method  (FEM) 
Analysis 


General 

An  elastic  response-spectrum  analysis  for  the  St.  Stephen  Powerhouse  was  performed 
to  evaluate  its  seismic  performance  using  the  displacement-based  procedure  discussed  in 
the  previous  chapters.  Development  of  the  discrete  finite  element  method  (FEM)  model 
for  elastic  response  analysis  of  the  erection  bay  is  discussed  below.  Two-dimensional 
(2D)  response  spectrum  analysis  was  performed  using  the  SAP2000  structural  analysis 
software.  The  erection  bay  was  analyzed  for  two  site-specific  motions  along  its  weakest 
direction,  namely  the  transverse  direction.  These  two  design  spectra  included  2,475-year 
equal  hazard  (probabilistic)  spectrum  and  the  MCE+a  determinisitic  spectrum.  Both 
spectra  are  described  in  Chapter  3.  Response  spectrum  analyses  were  conducted  to 
determine  the  effects  of  the  two  roof  support  boundary  conditions  described  earlier  and 
identified  as  “ideal”  “non-ideal”.  These  two  boundary  conditions  are  expected  to  envelop 
the  true  inertial  demands  on  the  powerhouse  superstructure  walls. 

FEM  Model 

A  discrete  FEM  model  of  the  generalized  transverse  section  of  the  erection  bay  is 
presented  in  Figure  4-1 .  The  base  of  the  substructure  is  assumed  fixed  at  EL  -41 .0.  The 
physical,  material,  and  boundary  conditions  for  the  model  are  described  in  Chapter  2. 

Through  various  parametric  studies  it  was  determined  that  the  behavior  of  the 
composite  substructure-superstructure  system  could  be  best  captured  using  beam-column 
elements  and  lumped  masses  to  represent  the  superstructure  and  shell  elements  and 
lumped  masses  to  represent  the  substructure.  The  beam-column  model  is  illustrated  in 
Figure  4-1. 

The  2-D  FEM  analyses  presented  herein  are  based  on  an  analytical  model  that 
represents  the  entire  erection  bay  in  terms  of  mass  and  stiffness.  The  substructure  was 
modeled  using  four-node  shell  elements.  In  general  the  substructure  is  the  region  below 
elevation  37.0  in  the  intake  area,  below  elevation  17.0  in  the  erection  bay  area  and  below 
elevation  57.0  in  the  tailrace  deck  area  (refer  to  Figure  1-3).  Superstructure  regions,  i.e. 
regions  above  the  substructure,  were  modeled  using  beam-column  elements.  Special 
beam-column  elements,  elements  that  could  approximate  the  stiffening  effects  of  the  two 
transverse  walls  located  in  the  intake  area  at  each  end  of  the  erection  bay,  were  used  in 
the  model.  These  special  elements  were  introduced  into  the  model  at  the  upstream  face 
between  elevation  37  and  elevation  86.  The  stiffening  effects  of  the  end  bay  wall  located 
at  the  south  end  of  the  erection  bay  was  not  included  in  the  model  since  its  influence 
would  be  limited  only  to  those  regions  of  the  erection  bay  located  in  close  proximity  to 
the  end  wall.  Masses  were  lumped  at  element  nodal  point  locations.  Upon  completion 
of  the  FEM  analysis,  the  total  earthquake  demands  on  various  wall  elements  were 
converted  to  force  demands  per  foot  of  wall.  This  was  accomplished  by  dividing  the  total 
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demand  by  the  erection  bay  width,  i.e.  66-feet,  This  action  was  taken  in  order  to  compare 
shear  and  moment  demands  with  shear  and  moment  capacities  determined  on  a  “per  foot” 
of  wall  basis. 


Figure  4-1 .  Beam-  column  model  of  the  erection  bay 

An  eigenvalue  extraction  of  the  lumped  column  model  indicated  that  a  few  of  the 
higher  modes  significantly  contribute  to  the  in-plane  motion  of  the  FEM  model, 
indicating  that  amplification  did  occur  in  the  superstructure  as  a  result  of  mass 
participation  by  the  substructure.  Modal  participation  factors  shown  in  Table  4-1 
indicates  the  significance  of  modes  6  and  8  in  amplifying  the  in-plane  motions  of  the 
erection  bay  model. 
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Table  4-1 


Modal  Participation  Factors  (MPF) 


Mode  No. 

Lumped  Model  | 

Period 

MPF 

1 

.301 

15 

2 

.2594 

15 

3 

.1072 

9 

4 

.080 

18 

5 

.068 

.8 

6 

.067 

14 

7 

.051 

9 

8 

.030 

33 

9 

.026 

4 

10 

.01 

2 

The  first  three  modes  of  vibration  for  the  erection  bay  along  with  the  higher  mode  8 
are  displayed  in  Figures  4-2  through  4-5.  Mode  8  illustrates  the  dominant  influence 
substructure  displacement  has  on  superstructure  response. 

Further  investigation  suggested  that  as  the  period  of  the  superstructure  wall 
approaches  the  period  of  the  substructure,  in  particular  for  wall  to  substructure  period 
ratios  between  one  and  two,  the  force  demands  on  the  wall  increase  significantly  over 
what  they  would  be  for  the  wall  acting  alone.  Higher  participation  factors  for  the 
substructure  modes  clearly  indicate  that  the  in-plane  flexibility  of  the  substructure  has 
significant  affects  on  the  superstructure  response. 
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Figure  4-2 


Erection  bay  first  mode  of  vibration 


Figure  4-3  Erection  bay  second  mode  of  vibration 
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Figure  4-4  Erection  bay  third  mode  of  vibration 


Figure  4-5  Erection  bay  eighth  mode  of  vibration 
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Earthquake  Demands 

The  elastic  response  spectrum  demands  for  “ideal”  and  “non-ideal”  support 
conditions  for  the  given  design  response  spectra  are  computed  using  the  lumped  mass 
beam-column  model  discussed  above.  The  assumptions  and  limitations  of  the  current 
analysis  are  discussed  below. 

•  A  lumped  beam  element  model  is  used  to  represent  substructure  and 
superstructure. 

•  Substructure  is  fixed  at  the  base 

•  Effective  widths  for  the  wall  elements  are  used  to  include  the  stiffness  reduction 
due  to  openings  on  walls. 

•  Elastic  response  analysis  using  SAP2000  was  based  on  LOAD-DEPENDENT 
RITZ  Method  (Acceleration  in  X-direction  only).  The  Ritz  method  is  proved  to 
be  a  faster  convergence  method  that  yields  more  accurate  results  than  the  use  of 
the  same  number  of  exact  mode  shapes  (Wilson,  1998).  The  RITZ  method 
computes  only  the  mode  shapes  that  are  excited  by  the  seismic  loading.  It  may 
require  including  over  40  modes  to  obtain  90%  participation  using  eigenvector 
method. 

•  The  out-of  plane  width  of  the  model  represents  the  entire  width  (66-feet)  of  the 
erection  bay 

•  The  effective  flexural  stiffness,  expressed  as  a  ratio  of  the  effective  moment  of 
inertia  to  the  gross  moment  of  inertia,  is  as  indicated  in  Table  1-1.  All 
substructure  elements  were  assumed  to  have  an  effective  stiffness  equal  to  80% 
of  gross  section  stiffness. 

•  The  damping  is  5%  and  the  design  spectra  are  the  2475-year  probabilistic  event 
and  the  deterministic  MCE  +  a  event. 

•  The  tee  beam  is  pinned  at  the  upstream  end  and  roller  at  the  downstream  end  for 
ideal  case  and  both  ends  pinned  for  non-ideal  case. 

•  The  crane  mass  and  load  is  distributed  equally  at  each  support. 

•  Dead  load  and  earthquake  loads,  i.e.  demands  from  the  response  spectrum 
analysis,  are  used  to  compute  the  total  demands  presented  in  Table  4-3. 

•  SAP2000  results  for  U  =  DL  EQ,  and  U  =  DL  -  EQ  are  identical.  The 
maximum  and  minimum  responses  for  both  cases,  however,  provide  the  positive 
or  negative  response  spectra  demands.  The  positive  or  negative  response  spectra 
demands  for  all  cases  are  listed  in  Table  4-3. 
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Figure  4-6  shows  the  critical  locations  on  the  erection  bay  model  for  which  the 
earthquake  demands  are  tabulated.  Modal  participation  factors  and  their  cumulative  sum 
for  both  ideal  and  non-ideal  cases  are  shown  in  Table  4-2.  This  table  indicates  that 
inclusion  of  first  10  modes  of  the  model  attains  over  95%  mass  participation.  Table  4-3 
shows  the  elastic  response  demands  for  critical  wall  locations  in  the  erection  bay.  In  this 
table,  the  maximum  and  minimum  responses  are  presented  for  a  load  combination  that 
includes  earthquake  loads,  dead  load,  and  live  loads,  i.e.  the  crane  weight  applied  equally 
as  a  distributed  load  on  both  upstream  and  downstream  walls  at  corbel  locations. 
Demands  include  the  moment,  shear,  and  axial  forces  for  the  entire  66-foot  wide  section 
of  the  erection  bay  model. 


Point  1 

Location  where  the  upstream 
wall  connects  to  the  intake  deck 
at  Elevation  86. 


1.5-Foot  Thick  Wall 


Point  3 

Location  where  the  downstream 
wall  connects  to  the  crane  rail 
corhel  at  Elevation  81.25. 


Point  2 


Location  where  the  downstream  wall  connects 
to  the  tailrace  deck  at  Elevation  57. 


Figure  4-6  Erection  bay  walls 
Critical  locations 
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Table  4-2 

MODAL  PARTICIPATION  FACTORS 


Mode 

Ideal 

Non-Ideal  | 

No. 

Period 

ux 

EUX 

Period 

UX 

lUX 

1 

0.301 

8.2 

8.2 

0.283 

16.4 

16.4 

2 

0.259 

8.6 

16.8 

0.115 

0.0 

16.4 

3 

0.107 

3.3 

20.3 

0.102 

7.1 

23.5 

4 

0.081 

12.3 

32.7 

0.075 

12.9 

36.1 

5 

0.070 

5.4 

38.1 

0.069 

1.6 

38.1 

6 

0.058 

2.5 

40.6 

0.057 

2.8 

40.9 

7 

0.049 

5.3 

45.9 

0.047 

5.9 

46.9 

8 

0.038 

6.4 

52.3 

0.035 

8.3 

55.1 

9 

0.027 

32.9 

85.2 

0.026 

30.9 

86.0 

10 

0.013 

10.4 

95.6 

0.012 

10.1 

96.1 

Table  4-3 

Member  demands  (SAP  2000)  for  both  NEHRP  2,474-Year  and  MCE  events  for 
“Ideal”  and  “Non-Ideal”  boundary  conditions  of  the  66  ft  wide  erection  bay 


Location 

Response 

Spectra 

Boundary 

Conditions 

Moment,  k-ft 

Max,  Min 

Shear,  k 

Max,  Min 

Axial,  k 

Max,  Min 

Point  1 

Elev.  86.00 

NEHRP  +DL+LL 

Ideal 

12,834 

(-12,579) 

1,604 

-116 

(-459) 

Non-Ideal 

18,716 

(-18,858) 

(2,397) 

-114 

J:461) _ 

MCE  +DL+LL 

Ideal 

7,498 

(-7,243) 

(931) 

-156 

(-419) 

Non-Ideal 

9,205 

(-9,347) 

(1,146) 

-155 

Point  3 

Elev.  81.25 

NEHRP  +DL+LL 

Ideal 

18,947 

(-19,202) 

(1,505) 

-242 

(-472) 

Non- Ideal 

7,928 

(-7,548) 

(636) 

-239 

(-475) 

MCE+DL+LL 

Ideal 

9,705 

(-9,960) 

(776) 

-251 

(-462) 

Non- Ideal 

4,115 

(-3,736) 

(335) 

-243 

Point  2 

Elev.  57.00 

NEHRP  +DL+LL 

Ideal 

48,174 

(-48,071) 

(2,100) 

135 

(-2,475) 

Non-Ideal 

25,294 

(-24,044) 

(1,614) 

-341 

(-1,968) 

MCE+DL+LL 

Ideal 

24,996 

(-24,893) 

(1,075) 

-434 

(-1,906) 

Non-Ideal 

13,503 

(-12,254) 

(850) 

-626 

(-1,684) 
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5  St.  Stephen  Erection  Bay  Walls 
Summary  of  Analysis  Results 


In  the  preceding  Chapters  the  process  for  obtaining  elastic  force  demands,  inelastic 
displacement  demands,  and  displacement  ductility  demands  based  on  the  results  of 
response  spectrum  analyses  was  described.  Methods  used  to  determine  force, 
displacement  and  displacement  ductility  capacities  were  also  described.  From  this 
information  demand  to  capacity  ratios  (OCR’s)  were  determined.  With  respect  to  shear 
and  sliding  shear  failure  mechanisms  demands  should  not  exceed  capacities,  or  the 
demand  to  capacity  ratio  should  be  less  than,  or  equal,  to  one  (OCR  <  1.0).  A  demand  to 
capacity  limit  of  one  is  placed  on  shear  and  sliding  shear  because  these  are  brittle  failure 
mechanisms.  A  potential  for  collapse  of  the  superstructure  exists  whenever  the  demand  to 
capacity  ratio  for  shear  and  sliding  shear  exceeds  one. 

The  demand  to  capacity  ratio  for  flexure  can  exceed  one  provided  the  displacement 
ductility  demand  does  not  exceed  the  displacement  ductility  capacity.  Or,  on  a  flexural 
displacement  ductility  basis  the  demand  to  capacity  ratio  should  be  less  than,  or  equal,  to 
one  (OCR  <  1 .0).  Whenever  the  flexural  displacement  ductility  demand  to  capacity 
exceeds  one,  rebar  splices  can  fail  or  reinforcing  steel  can  fracture.  Both  effects  can  lead 
to  collapse  of  the  superstructure.  The  locations  in  the  erection  bay  superstructure  where 
inelastic  action  will  take  place  during  major  earthquakes  as  described  previously  are 
designated  as  Points  1,  2,  and  3  and  illustrated  on  Figure  5-1 . 
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Figure  5-1  Erection  bay  superstructure 

Plastic  hinge  locations  (yield  regions) 


The  maximum  demands  on  these  plastic  hinge  regions  from  the  MDE  +  a,  2475-year, 
1000-year  and  475-year  earthquakes  are  presented  in  Tables  5-la,  5-lb,  5-lc,  and  5-ld 
respectively.  Demands  are  expressed  in  terms  of  moments,  shears,  and  displacement 
ductility.  Response  spectrum  analyses  for  the  1000-year  and  475-year  equal  hazard 
spectra  were  not  performed.  The  response  spectrum  analyses,  however,  are  linear  elastic 
and  the  equal  hazard  spectra  scale  uniformly  with  respect  to  all  periods  of  vibration. 
Therefore  the  demands  for  the  1000-year  and  475 -year  events,  for  a  given  period  of 
vibration,  can  be  estimated  based  on  the  ratio  of  the  spectral  acceleration  for  the  event 
under  consideration  (1000-year  or  475-year)  to  the  spectral  acceleration  for  the  2475-year 
event.  This  ratio  is  equal  to  0.45  for  the  1000-year  event  and  0.20  for  the  475-year  event. 
The  demands  for  each  location  (Points  1,  2,  or  3)  represent  the  maximum  demands  from 
the  “ideal”  or  “non-ideal”  boundary  condition  analytical  model.  The  non-ideal  condition 
produces  the  greatest  demands  at  Point  1 .  The  ideal  condition  produces  the  greatest 
demands  at  Points  2  and  3. 

The  capacities  of  the  various  members,  for  the  plastic  hinge  regions  illustrated  in 
Figure  5-1,  are  also  provided  in  Table  5-la  through  5-ld.  As  with  demand,  capacities  are 
expressed  in  terms  of  moments,  shears,  and  displacement  ductility. 

The  demand  to  capacity  ratios  (OCR’s)  for  each  plastic  hinge  region,  for  the  MDE  + 
(T,  2475-year,  1000-year  and  475-year  earthquakes  are  also  presented  in  Tables  5- la 
through  5-ld  respectively.  When  the  DCR  exceeds  the  allowable  value,  it  is  listed  in 
“bold”  text.  As  one  can  see  from  the  tables,  failure  mechanisms  that  can  lead  to  collapse 
of  the  erection  bay  superstructure  can  develop  for  all  but  the  475-year  event,  although  the 
number  of  potential  failure  mechanisms  decrease  as  earthquake  demands  are  reduced. 
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Force  demands  and  capacities  as  presented  in  the  tables  are  based  on  a  unit  length  of 
wall  (1-foot).  Moments  are  in  kip-feet  and  shears  are  in  kips. 

Although  the  displacement  demands  were  not  provided  in  the  tables,  the  inelastic 
displacements  at  the  roof  level  for  the  MCE  +  a,  and  2475-year  equal  hazard  event  were 
approximately  1-foot.  A  displacement  of  this  magnitude  will  place  severe  demands  on 
the  tee-beam  to  roof  connections  and  most  likely  lead  to  a  loss  of  support  for  the  tee- 
beams. 

The  maximum  lateral  force  demand  on  roof  fixed  bearings  is  approximately  equal  to 
the  maximum  shear  at  Point  1 .  From  the  shear  demand  information  in  the  tables  it  can  be 
seen,  that  except  for  the  475-year  event,  the  lateral  force  demand  will  exceed  the  fixed 
bearing  weld  capacity,  i.e.  1 1  kips  per  foot  (see  calculations  in  Chapter  1). 


Table  5.1a _ MCE  +  ct  earthquake 


Point  1  “Non-Ideal” 

Point  2  “Ideal” 

Point  3  “Ideal” 

Demand 

DCR 

Demand 

DCR 

Demand 

DCR 

9.9 

7.0 

1.4 

2.3 

3.0 

0.77 

4.6 

4.5 

1.1 

^■1 

Moment 

139 

33 

4.2 

379 

162 

151 

33 

4.6 

Shear 

17.4 

11.6 

1.5 

16.3 

27.0 

0.60 

11.8 

11.6 

1.0 

17.4 

17.5 

1.0 

16.3 

41.0 

0.40 

11.8 

17.5 

0.67 

Note:  Moment,  shear  and  sliding  shear  demands  are  per  foot  of  wall  and  are  equal  to  the 

demands  for  the  entire  erection  bay  (Table  4-3  values)  divided  by  the  erection  bay  width 
of  66-feet.  Moments  in  the  above  table  are  in  ft-kips  /  foot  of  wall.  Shear  and  sliding 
shear  are  in  kips  /  foot  of  wall. 


Table  5.1b _ NEHRP  2475-year  earthquake  demands 


Point  1  “Non-Ideal” 

Point  2  “Ideal” 

Point  3  “Ideal” 

Demand 

DCR 

Demand 

DCR 

Demand 

DCR 

liJHBliJW 

24.0 

7.0 

4.5 

3.0 

1.5 

8.8 

4.5 

2.0 

Moment 

286 

33 

8.7 

730 

162 

4.5 

291 

33 

8.8 

Shear 

36.0 

11.6 

3.1 

31.8 

27.0 

1.2 

22.8 

11.6 

2.0 

36.0 

17.5 

2.1 

31.8 

41,0 

0.78 

22.8 

17.5 

1.3 

Note:  Moment,  shear  and  sliding  shear  demands  are  per  foot  of  wall  and  are  equal  to  the 

demands  for  the  entire  erection  bay  (Table  4-3  values)  divided  by  the  erection  bay  width 
of  66-feet.  Moments  in  the  above  table  are  in  ft-kips  /  foot  of  wall.  Shear  and  sliding 
shear  are  in  kips  /  foot  of  wall. 
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Table  5.1c _ NEHRP  1000-year  earthquake  demands 


Point  1  “Non-Ideal” 

Point  2  “Ideal” 

Point  3  “Ideal” 

Demand 

DCR 

Demand 

DCR 

Demand 

DCR 

9.7 

7.0 

1.4 

2.0 

3.0 

0.67 

4.0 

4.5 

0.89 

Moment 

129 

33 

3.9 

329 

162 

2.0 

131 

33 

Shear 

16.2 

11.6 

14.3 

27.0 

0.53 

10.3 

11.6 

0.89 

Sliding  Shear 

16.2 

17.5 

0.93 

14.3 

41.0 

0.35 

10.3 

17.5 

0.59 

Note:  Moment,  shear  and  sliding  shear  demands  are  per  foot  of  wall  and  equal  to  the  values 
from  Table  5- 1  b  times  0.45.  Moments  in  the  above  table  are  in  ft-kips  /  foot  of  wall. 
Shear  and  sliding  shear  are  in  kips  /  foot  of  wall. 


Table  5. Id _ NEHRP  475-year  earthquake  demands 


Point  1  “Non-Ideal 

?? 

Point  2  “Ideal” 

Point  3  “Ideal” 

Demand 

DCR 

Demand 

DCR 

Demand 

DCR 

HIHH 

3.2 

7.0 

0.45 

0.9 

3.0 

0.30 

1.8 

4.5 

0.40 

Moment 

57 

33 

1.7 

146 

162 

0.90 

58 

33 

1.8 

Shear 

7.2 

11.6 

0.62 

6.4 

27.0 

0.24 

4.6 

11.6 

0.40 

7.2 

17.5 

0.41 

6.4 

41.0 

0.16 

4.6 

17.5 

0.26 

Note:  Moment,  shear  and  sliding  shear  demands  are  per  foot  of  wall  and  equal  to  the  values 
from  Table  5-lb  times  0.20.  Moments  in  the  above  table  are  in  ft-kips  /  foot  of  wall. 
Shear  and  sliding  shear  are  in  kips  /  foot  of  wall. 

It  must  be  recognized  that  the  shear  demands  indicated  in  the  above  tables  were 
obtained  from  the  elastic  response  spectrum  analysis  and  therefore  represent  the  shear 
demands  the  structure  would  experience  if  it  remained  elastic.  It  is  clear  from  the 
demand  to  capacity  ratios  for  moment  and  shear  in  the  above  tables  that  there  is  a 
hierarchy  in  the  formation  of  failure  mechanisms  and  that  flexural  yielding  will  occur 
before  shear  demands  reach  their  tabular  values.  In  reality  maximum  shear  demands  will 
be  limited  to  values  associated  with  the  onset  of  flexural  yielding.  The  shear  demands 
however  should  consider  the  potential  increase  in  moment  capacity  due  to  strain 
hardening  in  the  steel,  i.e.,  an  increase  in  steel  tensile  capacity  as  a  result  of  strain 
hardening  will  lead  to  an  increase  in  moment  capacity.  Normally  the  moment  capacity 
increase  due  to  strain  hardening  is  in  the  order  of  50-percent. 

Failure  will  likely  not  occur  due  to  shear  or  roof  bearing  connection  failure,  but  occur 
due  to  mechanisms  associated  inelastic  actions  in  plastic  hinge  regions.  This  means 
failure  will  occur  either  due  to  fracturing  of  the  flexural  reinforcing  steel  or  due  to  failure 
of  flexural  steel  rebar  splices  located  in  plastic  hinge  regions  (see  Table  2-1). 

Sample  calculations  illustrating  the  process  used  to  determine  displacement  ductility 
demand  and  displacement  ductility  capacity  are  illustrated  in  Table  5-2. 
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Table  5-2 

Sample  calculation  for  upstream  wall  displacement  ductility  capacity  Point  1  (EL  86) 


Displacement  Ductility  Demand  to  Capacity  Evaluation 
St  Stephen  PH  Erection  Bay 

Evaluate  the  1.5  foot  thick  section  of  the  upstream  wall 


April  25,  2000  File:  Ductility  1 .5  U/S  Foot  Wall 

Displacement  Ductility  Capacity 
(|)  := 0.0001 19  rads /inch 

^  From  M-(|)  Analysis 

(j)  y  1=0.002428  rads /inch 

Fy;=40.0 

d|j  1=0.875  inches 

Lp  ;=0.3Fy  dj, 

Lp=10.5  inches 
Lgff:=96.0  inches 


u^ 

|-  1  ♦ 

[Lp] 

1 

1 

p 

'"'’ll 

,l*yi 

\Leff/ 

L  ' 

ILeff/. 

p  =  7.019  Ductility  capacity 
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Table  5-2  (continued) 

Sample  calculation  for  upstream  wall  displacement  ductility  demand 
Point  1  (EL  86) 

2475-year  event  “non-ideal”  roof  bearing  support  condition 

NEHRP  TR  =  2475  Years 

Displacement  Ductility  Demand 


M  286 

Ft-Kips 

Mn-33 

Ft-Kips 

Tq:=0.20 

Seconds 

T§:=0.10  Seconds 

I^REQD'  — 

''^REQD= 

Priestley  Equation  for  Ductility  Demand 
n  dp'"^-^  (*^REQD“  ^)' 

^  DP  =  24  Use  for  ductility  demand 
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6  Conclusions 


A  performance-based  analysis  was  used  to  assess  the  seismic  performance  of 
St.  Stephen  powerhouse  erection  bay.  Two  site-specific  design  response  spectra  were 
used  to  evaluate  the  structure.  These  included  a  2,475-year  probabilistic  earthquake 
event  and  a  deterministic  MCE  +ct  event. 

The  St.  Stephen  Powerhouse  superstructure  is  vulnerable  to  collapse  for  earthquake 
hazards  with  return  periods  greater  than  500-years.  The  structure  does  not  meet  the  basic 
safety  objective  (collapse  prevention)  of  FEMA  356. 
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MEMO  FOR  CECW-E 
CECW-0 

SUBJECT:  Seismic  Stability  of  St.  Steven's  Hydropower  Plant,  South  Carolina 

1.  Summary.  At  the  request  of  CESAD  and  CESAC,  a  meeting  was  held  to  review  the 
seismic  stability  of  St.  Steven's  Hydropower  Plant,  South  Carolina.  In  conjunction  with  a 
presentation  on  the  development  of  design  ground  motions  and  structural  analysis  of  the 
project,  it  was  determined  that,  based  on  existing  guidance  in  ER  1 110-2-1806,  the  St. 
Steven's  Hydropower  Plant  is  not  a  "critical"  structure. 

2.  Place  and  Date.  Information  Technology  Laboratory,  Waterways  Experiment  Station, 
Vicksburg,  Mississippi,  2  May  2000. 

3.  Attendees:  Jerry  Foster,  CECW-ET,  Michael  Klosterman,  CECW-EG;  Byron  Foster, 
CESAD-ET-E,  Wayne  Bieganowsky,  CESAC-TS-DF,  Robert  Ebeling,  CEERD-ID-P, 
Wayne  Jones,  CEERD-ID-P,  Robert  Hall,  CEERD-SS-A,  Mostafiz  Chowdhury,  CEERD- 
SS-A,  Don  Yule,  CEERD-GG-H,  Dr.  Collins,  FERC  (conference  phone),  Ralph  Strom, 
consultant,  and  Wayne  Jones,  CEERD-ID-P,  visitor. 

4.  Purpose.  To  review  the  seismic  stability  of  St.  Steven's  Hydropower  Plant,  South 
Carolina. 

5.  Background.  St.  Steven's  Hydropower  Plant  is  located  on  the  Santee  River  and  is  part 
of  the  Pinopolous  Reservoir  system  near  Charleston,  South  Carolina.  The  project  was 
completed  in  1985  and  is  being  reevaluated  for  seismic  stability  on  a  15-year  cycle  in 
accordance  with  ER  1 110-2-1806,  Earthquake  Design  and  Evaluation  for  Civil  works 
Projects.  Significant  advances  in  the  state-of-the-art  of  earthquake  engineering  have 
occurred  since  the  project  was  completed  and  led  the  district  to  believe  that  the  dam 
would  be  a  candidate  for  remediation  under  the  Dam  Safety  Assurance  Program.  Design 
ground  motions  and  structural  analysis  of  the  project  were  performed  based  on  these 
recent  advances.  Interpretation  of  ER  1 110-2-1806  led  the  district  and  CEERD  personnel 
to  conclude  that  the  St.  Steven's  Hydropower  Plant  was  a  critical  structure  and  subject  to 
the  new  NEHRP  guidelines  (FEMA-273  &  FEMA-274). 

6.  Discussion. 

a.  Criticality  of  the  structure. 

(1).  According  to  guidance  contained  in  ER  1 1 10-2-1806,  para.  5.a.  and 
Appendix  B,  "Critical  features  are  the  engineering  structures  ...whose  failure 
during  or  immediately  following  an  earthquake  could  result  in  loss  of  life.  Loss  of 
life  potential  (is)  based  upon  inundation  mapping  of  (the)  area  downstream  of  the 
project."  The  most  likely  failure  scenarios  developed  for  the  St.  Steven's 
hydropower  plant  indicate  there  is  no  incremental  loss  of  life  downstream  of  the 
structure  due  to  its  failure  during  an  earthquake  event. 
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(2) .  Other  dams  in  the  Santee-Pinopolous  Reservoir  system  control  the 
consequences  of  failure  that  would  occur  during  a  major  seismic  event.  The  failure  of 
St.  Steven's  powerhouse  would  have  no  incremental  impact  on  loss  of  life  downstream. 

(3) .  ER  1 110-2-1806,  para.  8.b.  indicates  that  "Minimum  standards  for  the  seismic 
design  or  evaluation  of  (occupied)  buildings  ...are  available  in  national,  regional,  or  local 
building  codes."  The  presence  of  USACE  personnel  in  a  structure  does  not  qualify  the 
structure  for  the  'critical'  designation  of  para.  5.a.  of  ER  1 110-2-1806.  The  National 
Earthquake  Hazard  Reduction  Program  (NEHRP)  is  designed  to  insure  the  safety  of 
personnel  in  USACE  and  other  Federal  buildings  in  accordance  with  PL  101-610  and 
Executive  Order  12941.  The  seismic  safety  of  all  USACE  buildings  has  been  evaluated 
under  this  program.  The  district  should  verify  the  seismic  risk  assigned  to  the  St.  Steven's 
hydropower  plant  during  these  evaluations.  The  program  proponent  is  Helen  Petersen, 
CENWD  (503-808-3833). 

b.  Appropriate  ground  motions.  Probabilistic  ground  motions  are  obtained  from 
ER  1 1 10-2-1806,  Appendix  D.  Appendix  D  consists  of  1991  USGS  seismic  maps;  a  1998 
revision  has  updated  this  appendix  to  the  1995  USGS/NEHRP  seismic  maps.  The  most 
recent  1997  NEHRP  seismic  maps  have  not  been  adopted  for  use  in  ER  1 1 10-2-  1806.  At 
the  present  time  ground  motions  used  to  evaluate  civil  works  projects  will  be  obtained 
using  the  1995  USGS/NEHRP  seismic  maps. 

c.  Future  revisions  to  ER  1 1 1  0-2-1806  may  incorporate  concepts  addressed  in 
FEMA  273  and  274.  At  the  present  time,  all  policy  and  guidance  regarding  earthquake 
design  and  evaluation  for  civil  works  projects  is  contained  in  the  existing  ER  1 110-2- 
1806. 


7.  Future  Action.  The  St.  Steven's  Hydropower  Plant  should  be  evaluated  as  a 
noncritical  structure  according  to  the  requirements  of  (ER  1 1 10-2-1806.  It  should  be 
demonstrated  that  the  structure  can  withstand  the  MDE  for  noncritical  structures  (ER 
1 1 10-2-1806,  para.  5d  and  para.  6b.  Any  justification  for  proposed  remediation  to 
withstand  seismic  loads  in  excess  of  the  MDE  must  be  based  upon  an  economic  analysis 
in-accordance  with  guidance  for  the  Major  Rehabilitation  Program. 


JERRY  FOSTER,  PE  MICHAEL  J.  KLOSTERMAN,  RPG 

Structural  Engineer,  CEC W-ET  Chief  Geologist,  CECW-EG 
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MEMORANDUM  FOR  Dr.  Robert  Ebeling,  CEERD-ID-P,  US  Army  Engineer  Research  and 
Development  Center,  3909  Halls  Ferry  Road,  Vicksburg,  MS,  39180-6199 

SUBJECT:  Corps  of  Engineers  Hazard  Reduction  Program  -  Powerhouse  Evaluation  Criteria 

1 .  The  Corps  of  Engineers  Hazard  Reduction  Program  (CEHRP)  was  implemented  to  assure  the 
safety  of  all  occupants  in  USACE  owned  structured  in  accordance  with  Executive  Order  12941 
and  Public  Law  101-614,  National  Earthquake  Hazard  Reduction  Program  (NEHRP).  As 
required  in  EO  12941,  the  Interagency  Committee  on  Seismic  Safety  in  Construction  (ICSSC) 
adopted  minimum  life  safety  standards  for  all  federally  owned  buildings  and  building  type 
structures  based  on  existing  Federal  Emergency  Management  Agency  (FEMA)  guidance 
documents. 

2.  In  complying  with  the  minimum  standards  adopted  by  the  ICSSC  and  the  intent  of  those 
standards  the  drawings  of  many  of  the  USACE  owned  powerhouses  were  reviewed,  the  districts 
responsible  for  those  powerhouses  were  surveyed,  a  number  of  studies  were  performed,  and 
existing  seismic  criteria  was  reviewed.  With  few  exceptions,  powerhouse  superstructures  were 
considered  non-critical  features.  Specifically,  failure  of  the  powerhouse  superstructure  is  not 
likely  to  result  in  failure  of  the  dam  and  consequently  loss  of  life  downstream.  This  means  that 
site  specific  earthquake  assessments  are  not  required  and  are  only  recommended  if  it  is  likely  that 
the  site  specific  assessment  will  reduce  the  earthquake  demand  on  a  structure.  In  addition,  it 
allows  the  use  of  progressively  more  complex  evaluation  procedures  as  the  specific  situation 
warrants. 

3.  However,  for  those  powerhouses  not  remotely  operated,  no  reasonable  justification  to 
establish  a  performance  level  or  evaluation  criteria  less  than  that  required  for  any  occupied 
building  could  be  found.  To  comply  with  the  intent  of  the  minimum  ICSSC  standards,  the 
FEMA  building  evaluation  criteria  adopted  as  the  minimum  life  safety  standard  (FEMA  178, 
FEMA  310)  were  reviewed  and  the  recommended  rehabilitation  criteria  (FEMA  273,  FEMA 
274,  and  FEMA  356)  were  consulted  and  consolidated  to  a  single  draft  criteria  document  with 
progressively  more  complex  evaluation  procedures. 

4.  The  powerhouse  superstructure  evaluations  conducted  and  fimded  by  our  office  were  based  on 
a  design  level  earthquake  with  a  collapse  prevention  check  for  the  maximum  considered 
earthquake.  The  design  level  earthquake  used  was  approximately  a  1000  year  probabilistic 
earthquake  event  based  on  the  1997  NEHRP  maps.  In  addition,  in  accordance  with  FEMA  273 
and  FEMA  274  a  collapse  prevention  check  was  performed  using  the  Maximum  Considered 
Earthquake  based  on  the  1997  NEHRP  maps.  On  these  maps  the  Maximum  Considered 
Earthquake  corresponds  to  either  a  2475-year  event,  or  in  zones  near  major  active  faults  (i.e. 
California)  to  either  the  smaller  of  the  2475-year  event  or  150-percent  of  the  mean  ground 
motion  for  a  deterministic  event  (i.e.  Maximum  Credible  Earthquake).  The  decision  to  use  of  the 
1997  NEHRP  maps  instead  of  the  maps  shown  in  ER  1 1 10-2-1806  was  based  on  the  fact  that  this 
was  the  latest  balloted  consensus  version,  the  recommendations  of  several  industry  consultants 
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some  of  whom  were  involved  in  the  development  of  the  maps,  and  the  FEMA  reporting 
requirements.  The  seismicity  of  all  buildings  and  structures  in  FEMA  360  (the  report  to 
Congress)  were  required  to  reflect  the  1997  NEHRP  maps. 

5.  The  evaluation  of  St.  Steven’s  Hydropower  Plant  in  South  Carolina  performed  by  our  office 
for  the  design  level  earthquake  described  above  resulted  in  the  conclusion  that  the  powerhouse 
does  not  meet  the  recommended  performance  level,  the  minimum  life  safety  standard  required  by 
the  ICSSC.,  or  the  basic  safety  objective  outlined  in  FEMA  356  and  TI  5-809. 

6.  Based  on  discussions  with  Dr.  Ebeling  and  his  review  of  the  criteria  used  in  our  evaluations, 
the  seismic  stability  evaluation  conducted  by  CEERD  of  the  St.  Steven’s  Hydropower  Plant  is 
consistent  with  the  FEMA  356  (basic  safety  objective,  collapse  prevention)  seismic  analysis 
procedures  and  consistent  with  the  methodology  used  in  the  CEHRP  powerhouse  evaluations. 


Ms  Helen  Petersen 
Structural  Engineer 
CEHRP  program  Manager 
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